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Semiconductor p + -p~-n-p + -n ++ structures with large junction and con- 
tact areas are treated as 1 x 2-dimensional active media, in which self-organized 
pattern formation can be expected. The local bistable behavior of the struc- 
tures may emanate from two different mechanisms both governed by a nonlin- 
ear current feedback-loop between the electrons and holes injected from the 
outer layers. By considering the device to be composed of an active subsystem 
with negative differential resistance and a passive resistive layer with positive 
differential resistance an analytical approach is suggested to understand and 
describe the corresponding physical mechanisms in a self-consistent way. Ana- 
lytical solutions of the derived model equations allow a description of homoge- 
neous stationary states and yield explicit expressions of the current-density vs. 
voltage characteristics of the whole structure and its subsystems. A stability 
analysis of the homogeneous states with respect to two-dimensional transver- 
sal harmonic fluctuations is performed for the two cases under study. The 
resulting dispersion relations allow two different types of instability. While 
the first one is of Ridley's type which is characteristic for any spatially ex- 
tended electrical system with negative differential resistance, the second type 
can be considered as a solid-state analogue of Turing's instability known as 
a generic instability mechanism which may lead, e.g., to the formation of 
periodic patterns. 



I. INTRODUCTION 

Basic ideas on instabilities and dissipative pattern formation in open spatially extended 
nonlinear systemsllhi have gained substantial interest in many fields of science in the last 
two decades, (see, e.g., Refs-BifliU) . The investigations become rather intensive also in the 
field of solid-state physics and electronicsBi'0 where a large variety of different nonlinear 
nonequilibrium phenomena occurs, which can cause spontaneously arising spatio-temporal 
patterns in materials and devices. Together with highly developed device technology, these 
circumstances suggest promising opportunities for future research and developements. 

Recently, attention has been called to thyristor-like semiconductor structures with large 
active areasEJllil, as these nonlinear systems with bistable properties show several spatial and 
spatio-temporal current-density patterns. Such semiconductor structures could potentially 



be used as multi-stable elements for integrated circuits, self-organizing devices for image 
recognition and image processing, etc.iE10. One of the most striking results of former con- 
siderations based on phenomenological mo delHl is that the instability mechanisms posess 
some features very similar to those studied in biology or chemical media, which may show an 
extremly manifold selforganizing behavior. However, at present the understanding of insta- 
bility mechanisms and pattern formation in multilayer structures is far from being complete. 
Apparently, there is a lack of systematic investigations concerning self-organization aspects 
for this family of devices. In a very early theoretical paper0 the process of single current 
filament formation has only been analyzed for the simplest case of p-n-p-n structures, that 
according to Ridley's terminologyEl belong to the class of extended electrical systems with 
current controlled negative differential resistance. Recently, gate-driven p-n-p-n structures 
have been used to study current filamentation and the dynamics of switching fronts under the 
constraints of two global couplingsEli^. In specially designed multilayer structures different 
types of filament oscillation have been observedilHil. Very few theoretical papers have been 
devoted to instabilities in standard thyristors; they all aimed merely on a safe operation in 
power applications^!^. Furthermore, it is not an easy task to adapt the advanced experi- 
ence of power device technology directly to the specifics of devices showing phenomena of 
self-organization. One of the main reasons for this is that the operation of semiconductor 
devices and especially thyristors depends very sensitively on many physical and design pa- 
rameters. On the one hand, this gives numerous opportunities for a proper choice of those 
parameters responsible for the intended function of these devices. But on the other hand the 
high-dimensional parameter space makes a systematic analysis rather difficult, because the 
use of numerical methods and the performance of experiments are usually time-consuming 
and/or expensive. Therefore detailed studies of the underlying physical processes in con- 
junction with the material and the design parameters are necessary and should lead to the 
deepest possible analytical description of the device physics. 

Following this motivation and aiming at further development of concepts introduced ear- 
lier in Refs. we suggest in this paper an analytical approach for instability mechanisms 
in a family of devices based on thyristor-like p + -p~-n-p + -n ++ structures, in which several 
modes of pattern formation have already been found experiment ally0. 

The outline of the paper is as follows: In section 2 the investigated semiconductor 
structure and its basic physics with emphasis on the nonlinear properties are described. 
A non-stationary 1 x 2-dimensional model rendering possible an analytical investigation of 
two-dimensional current-density and potential distributions transversal to the main current 
flow direction is derived in section 3 for two cases characterized by different nonlinear elec- 
trical mechanisms. In section 4, stationary homogeneous states are determined from the 
model equations and used to derive the current- voltage characteristics for systems with uni- 
form current-density distributions. Finally, the stability of stationary homogeneous states is 
studied in section 5, including a discussion of important parameter dependencies. 



II. MULTILAYER SEMICONDUCTOR STRUCTURES AS lx 2-DIMENSIONAL 

EXTENDED ACTIVE MEDIA 

The system under consideration consists of a semiconductor structure with a vertical 
width w of several hundreds of /zm and with the two other (transversal) widths L y ,L z ^> w. 
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Along the vertical direction the structure is characterized by a p + -p~-n-p + -n ++ doping profile 
as schematically shown in Fig. [l|. The outer p + and n ++ layer of the structure are provided 
with metal contacts. The internal design of the structure possesses the following peculiarities: 

(i) Along the (vertical) x-direction the p + -p~-n-p + -n ++ structure can be divided into an 
active (triggering) subsystem consisting of a four-layer thyristor-like p~-n-p + -n ++ structure 
and a subsystem of a distributed resistive layer consisting essentially of the remainder p~ 
layer. Note that the thyristor-like regeneration mechanism in the active subsystem causes 
an autocatalytic increase of the current density in a certain current interval which is coun- 
teracted by the resistive layer. 

The width of the n-p + -n ++ part is only a few tens of /mi which is more than an order 
less than the whole width w of the device. The lifetime r of the excess carriers in the bulk 
material is supposed to be so small that conductivity modulation of the bulk due to the 
injection of electrons and holes from the outer n ++ and p + layers is prevented. This can be 
achieved by reducing the lifetime in the main part of the bulk in a certain distance to the 
anode contact leaving a relatively thin part of the p~ layer with large r nearby the inner n 
layer. 

(ii) The doping concentrations Nj[ and Np in the layers forming the cathode emitter p + - 
n ++ junction [Fig. [TJ(b)] are similar to those used in high frequency bipolar transistors. The 
electron injection efficiency of this junction increases from nearly zero at very low current 
densities to a saturation value of approximately 0.8-0.9 at rather small current densities of 
the order of 1 mA/cm 2 . 

(iii) The doping concentration of the p~ layer is much smaller than the doping 
concentration Np of the n base. Hence, for low currents, current transport through this 
junction is dominated by a leakage current of electrons from the n base into the p~ emitter 
and, consequently, the injection efficiency of the p~-n junction is very small. 

(iv) The widths of the p + and n base, w\ and W20, in the thyristor-like subsystem are 
much smaller than the diffusion lengths L n and L p of the minority carriers in these layers: 
W\ <C L n , W20 <C L p . Consequently, both base transport factors are about one. 

When a dc voltage is applied to the device such that the outer n ++ layer is biased 
negatively with respect to the outer p + layer, the anode and cathode emitter junction are 
forward biased, while the collector junction is reverse-biased. Thus, double injection occurs 
inside the thyristor-like subsystem accompanied by a characteristic nonequilibrium plasma- 
field stratification along the vertical direction x as shown in Fig. 0(a). Due to the chosen 
doping proportion at the p~-n junction, which is just inverse to that in standard thyristors, 
there is a significant electron leakage current from the n base into the p~ layer which obviously 
exceeds the electron injection current. Therefore, a rather large excess charge accumulates 
in a thin plasma layer (V layer). If the carrier lifetime in this plasma layer is sufficiently 
large, the local dynamics of this charge at non-stationary conditions may be considered to 
be the most inertial process, thus determining the temporal evolution of the whole device. 
The derivation of the model equations in section 3 is based essentially on these physical 
properties. 

The nonlinear feedback mechanism between the electron and hole injection currents in the 
p~-n-p + -n ++ subsystem is similar to that in standard thyristors. However, there are remark- 
able peculiarities and it is useful to distinguish between two cases, corresponding to certain 
current intervals, which are characterized by a particular nonlinear mechanism, respectively, 
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causing the appearance of negative differential resistance in the current- volt age characteris- 
tic of the thyristor-like subsystem. In analogy to the two-transistor approach for thyristors 
we consider the thyristor-like subsystem to be composed of an n-p + -n ++ and a p~-rwp + 
transistor to specify the two different regenerative mechanisms. Current regeneratiorHEa in 
thyristors is possible for current intervals, in which the condition = dVr/dj < holds, 
where Vr and j denote the thyristor voltage and current, and is the differential resistance. 
dVr/dj is implicitly given by the steady state condition 

ai (j,V T ) + a 2 (j,V T ) = 1, (1) 

where a\ and a 2 denote the current gains of the n-p + -n ++ and the p~-n-p + transistor, 
respectively [Fig. 0(a)]. 

Charge carrier recombination in the space charge region of the emitter junction of a 
transistor influences the current gain significant lycalH3 and typically leads to a monotonic 
increase of the current gain with increasing current for not too large currents. 

Consider the device parameters to be of such kind that the current gain ot\ approaches 
a value close to one at relatively low currents (j < 10 _3 A/cm 2 ). Then Eq. (P can be 
fulfilled at a very low injection efficiency of the p~-n emitter junction. Thus, this case - 
subsequently called low current range - is roughly characterized by < «i = ai(j) < 1 and 
< a 2 ~ const 1. 

If the doping profile and other parameters are chosen such that a± saturates in the low 
current range before Eq. ([!]) is true, regeneration does not occur at low currents. Nev- 
ertheless, at larger currents regeneration might become possible again when a 2 increases 
monotonously with the current. This can occur at sufficiently large current, when a plasma 
layer V evolves in the p~ layer near the p~-n junction and the injection of excess holes from 
this plasma layer into the n base increases superlinearly as proposed by FletcheiR Conse- 
quently, the V layer plays the role of a plasma emitter, the efficiency of which rises strongly 
with the current so that Eq. ([I]) can be fulfilled. This case is called moderate current range 
subsequently and is roughly characterized by < «i ~ const < 1 and < a 2 = a 2 (j) < 1. 

Before we start with a detailed analysis of these two cases let us turn to a short discussion 
on important transversal processes which have to be taken into account for an analytic 
description, too. To this end, consider a so-called differential element of the whole p + -p~-n- 
p + -n ++ structure, which is small in y and z direction. Fig. 0(b) shows schematically contours 
of possible current fluctuations in such an element which may lead to a destabilization of 
the homogeneous current flow. Obviously, there are three channels along which transversal 
currents effect a transversal spatial coupling in the thyristor-like subsystem: the p + base, 
the n base, and the possibly conductivity modulated V layer in the p~ bulk, in which 
the corresponding transversal fluctuations 8i\, 5i 2 , and Si-p may arise. For the resistor 
subsystem transversal current fluctuations 51r are taken into account; they are distributed 
over the whole bulk due to three-dimensional deformations of the potential distribution 
(p(x, y, z, t); these deformations are caused by potential fluctuations at the interface between 
the two subsystems, which in turn result from local current fluctuations in the thyristor-like 
subsystem. 

The metallic layers forming the anode and cathode contact on the top and bottom of the 
device represent further channels denoted by 5ia and Six in Fig. 0(b) and render possible 
the formation of closed loops of current fluctuations inside the system without any inter- 
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action with elements of the external circuit. Thus, all differential elements of the system 
are connected to each other by both the vertical feedback mechanism and transversal cou- 
plings which have to be considered in a self-consistent description as outlined in detail in the 
following sections. 

III. NON-STATIONARY (lx2)-DIMENSIONAL ANALYTICAL MODEL 

A. The thyristor-like subsystem 

As mentioned above the excess plasma of the p~ layer is considered to be the most inertial 
part of the system. Consequently, all other nonstationary transport effects are neglected. 
We start with a derivation of the basic equations for the low current range. Based on these 
results the peculiarities concerning the case of moderate currents are discussed in a second 
subsection. 

1. Low current range 

In the low current case the dynamics of the plasma emitter can be described in terms of 
the excess charge density Q = Q(y, z, t) per unit square in the plasma layer V: 

dQ Q 

T7~ = h jln ~ 3Tn + V_l • t Vn , (2) 

at r 

where 

ivn = ~b{< n > /N A )a P w v V ± V T + D n V ± Q « D n V±Q (3) 

is the sheet electron current-density in the plasma layer, Vj_ = (d/dy,d/dz) the two- 
dimensional V-operator, D n the diffusion coefficient of electrons, r the lifetime of excess 
charge carriers at low injection, j 2n the x component of the electron current-density at the 
boundary between the plasma layer V and the n base, jxn the x component of the electron 
current- density leaving the thyristor-like subsystem and entering the resistor subsystem.^] Vr 
denotes the voltage drop across the thyristor-like subsystem, N A the acceptor concentration 
of the p~ layer, a> the conductivity of the plasma layer V, which is equal to qfi p N A in the 
low current range, q the elementary charge, b = fi n /fi p , fi n and fi p the electron and hole 

W-p 

mobility, and < n >= j n(x)dx the average electron density per unit square in the plasma 

o 

layer with width w-p. 



1 Note that all quantities denoted by j... refer to the x component of the corresponding current- 
density vector; for a better readability the supplement x component is omitted in the following 
explanation. Accordingly, quantities denoted by i... = (i..., y ,i..., z ) refer to pure sheet current den- 
sities per unit length in a certain layer of the structure. Physical constants as well as material and 
design parameters are listed separately in Tables | and O. 
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Following the charge-control model of transistors and thyristors (see, e.g., Ref.t3) the 
variable jx n can be expressed by 

Q n W V fA\ 



' n 



Thus, the charge equation can be re-written as 

?Q- = D n A ± Q + j 2n -Q-, (5) 
at r* 

with 1/t* = 1/t+1/9. The variable Q can be further approximated in the following way by 
means of the maximum value of the excess charge carrier concentration in the plasma 
layer: 

Q K TO?*. (6) 

Relations connecting the concentrations of the minority charge carriers at both sides of the 
jT-n junction, rift and p™, and the hole current- density j 2p at the interface between the 
plasma layer and the n base are provided by the classical theory of p-n junctions (see, e.g., 
RefS): 

r$NZ=i%N D , (7) 

<l D pPn / Q x 



JCp = J2p 



w 2 {V c 



c 



where is the donor concentration in the n base, w 2 the effective width of the n base, 
which depends on the voltage drop Vc across the n-p + collector junction, and D p the diffusion 
coefficient of holes. Note that due to the assumption w 2 <^ L p [Sect. 2, assumption (iv)] the 
hole current- density jc P at the collector is equal to j 2p . 

The dependence of thebase width w 2 on the collector voltage drop Vc(y,z,t) can be 
approximated by (see Ref.t^l): 



w 2 (V c ) = w 2Q 



1 - 



■,l rr qN D W 2Q 

with V p th = — 

2€£n 




(9) 



where w 2 o denotes the whole width of the n base, V pt h the punch-through voltage of the n-p + 
junction, eo and e the vacuum permittivity and the dielectric constant of silicon. 

Equations (^|) - (H) yield an expression linking the hole current-density at the collector 
junction, j Cp = j 2p , the charge Q and the collector voltage V c : 

(10) 

\ V v pth / 

AT- 9 n 

with , . , 

N D wpw 20 
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In the n-p + -n ++ transistor the local current balance for a differential element of the p + 
base is given by: 

jc - ji = V_l • ii = -erii/;iA.i.</?i, (11) 
with (Ji = qHpN^, 

where jc and jx denote the total current-densities at the n-p + collector and p + -n ++ emitter 
junction, o\ the conductivity of the p + base, and (fx the voltage drop across the emitter 
junction. 

Because of u>i C L„, the electron current-density at the collector j Cn = is practically 
equal to the electron current- density j\ n injected by the n ++ layer. Nevertheless, since the 
hole leakage from the p + base into the n ++ emitter is strongly nonlinear we have to take 
into account a j-dependence (or bias-dependence) of the parameter a±. The parametric 
dependence between the total emitter current-density jx and its electron injection component 
ji n is given as an implicit function coupling the emitter voltage tp% with the emitter current- 
density jx'. 

= F(jx,ip) = (jxns + ji P s)exp (qipx/kT) + j 1Rs exp (Aipx) - jx, (12) 
jcn ~ jin = jins exp (qipx/kT). (13) 

Here jx ns , jips and jxR S are saturation values of the electron injection current- density, and the 
linear and nonlinear hole leakage current-densities, respectively. The .parameter A depends 
on the concrete mechanism of the hole leakage in the p + -n ++ junctions. For the case under 
study we assume A = q/2kT, i.e. Sah-Noyce-Shockley recombination in the space charge 
region. 

In the V-n-p + transistor the local current balance for a differential element of the n base 
and of the plasma layer V is given by: 

h ~ Jc = V_l • i 2 

= J2n — jCn 

= -V ± [a 2 w 2 (Vc)V ± (V c + (fx)) , (14) 
3t ~ 32 = V_l • Vp 

= -a v w v A ± V T + D n A ± Q, (15) 

(16) 

where a 2 = qfi n N D is the conductivity of the n base. Note, that low injection conditions are 
assumed for the plasma emitter. 

Combining Eqs. ( |14|) and dig) yields: 

jc = 3t - V_l • i-p - V_l • i 2 

= 3t + [a 2 w 2 (V c )V j_(V c + fx)} + a P w v A ± V T - D n A ± Q, (17) 
jc P = jc ~ jcn- (18) 

The voltage drop Vr across the thyristor-like subsystem is given by Vr = Vq + <fi + <P2, 
where <p 2 is the voltage drop across the p~-n junction. Using p™ = p n0 exp (q<p 2 /kT) = 
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n p o(N A /NE>) exp (qip 2 /kT) and Eqs. @ and (0) the following relation between Q and (^2 can 
be derived: 

Q = exp (qip 2 /kT) 

= ^^exp( q( p 2 /kT), (19) 

where p n o and n p o denote the equilibrium minority carrier concentrations of the n base and 
p~ bulk, respectively, and rii is the intrinsic carrier concentration. 



2. Moderate current range 

Now we indicate the peculiarities occuring at moderate currents, which case is charac- 
terized by a high injection level in the p~ bulk. In this case ambipolar effects in the plasma 
layer should be included into consideration. In the frame of the charge-control model the 
following relations between the current densities jTn and jr, the sheet current densities i-p n 
and i-p = i-p n + *-pp an d the charge Q in the plasma layer V are valid: 

b Q 

= (T+T + s- < 20 > 

i P = i PTl + i Vp = -w v a P (Q)V ± V T + (b - 1)D P V ± Q, (21) 
im = —^w v cr v (Q)V ± V T + bD p V ± Q. (22) 

Here w-p<r-p(Q) = (6 + l)p. p Q is the transversal conductance of the plasma layer, 6^ = 
w v /(2Dh) the transit time through the plasma layer, and Dh = 2bD p /(b+ 1) the ambipolar 
diffusion coefficient. Note that Eqs. ( pOl) and (p2|) substitute Eqs. @ and (H), respectively. 
The charge balance equation for the plasma layer should be re-written as 

-WT = h J2n - jTn + Vj_ • l<p n , (23) 

Ot T h 

where is the high-injection lifetime. By inserting Eqs. (|20|) and (p2|) we obtain 

^ = o£ p A ± Q - 9- - 9. + j 2n - JL^j T - ± Vl [w v a v (Q)V x V T ] (24) 
at Th Uh o+i o+i 

[compare with Eq. (|)]. 

The relation between the excess carrier concentrations at the borders of the space charge 
region of the p~-n junction changes from Shockley's to Fletcher's^ form: 

(n£) 2 = p™N D (25) 

[compare with Eq. ([/])]. 

Taking into account Eqs. (||), (||), and (^5]) we obtain in analogy to Eq. ( PH) an expression 
linking the charge Q in the plasma layer, the voltage drop Vr ~ Vc across the V-n-p + -n ++ 
subsystem and the hole current- density jc P of the collector: 



8 



w 2 {V c ) 


1/2 




1/2 


w 2 o 




(1 - ai)j m 





Q = Qr 



where the normalizing quantities Q m and j m are connected by 



(26) 



Qr. 



qN D w^W2o(l - ai)j r , 



.1/2 



(27) 



As already mentioned above, the current gain a\ of the n-p + -n ++ transistor is approx- 
imately current-independent in the moderate current range. Then, the electron current- 
density jcn of the collector is given by 



JCr, 



with oi\ = const > 



1 



(28) 



In the V-n-p + transistor conductivity modulation in the plasma layer V due to excess 
carrier concentrations should be taken into account at high injection level; consequently 
Eq. (|15) is substituted by: 

3t ~ 32 



Vj_ • i-p 

(b - 1)D P A L Q 



V ± [w v a v (Q)V ± V T ] . 



(29) 



On account of vanishing field and diffusion currents at the transversal borders of the sam- 
ple (y = 0, L y ; z — 0, L z ) the following boundary conditions for the V-n-p + -n ++ subsystem 
have been chosen for both the low-injection and the high-injection case: 



(TO • V±Q)\y= ,L y ; z=0,L z = 0, 
(TO • V x<^l)\y=0,L y ; z=0,L 2 = 0, 
(TO • V_LV? 2 )|y=0,L H ; z=0,L z = 0, 
(TO • Vj.^)| tf== o,L wi z=0,L z = 0. 



(30) 
(31) 
(32) 
(33) 



Here n is a unit vector normal to the transversal surfaces of the sample. 



B. Resistor-like subsystem and external circuit 

In the resistor-like part of the p~ bulk displacement currents are neglected. The three- 
dimensional potential distribution in this region is described by Laplace's equation: 

Atp = (34) 

with the boundary conditions 

<p\ x =o = Va, <p\x=w r = V T , (35) 

(TO • V ±(p)\ y=0 ,L v ; z=0,L z = 0, (36) 

where Va denotes the voltage drop across the whole sample, wr the width of the resistive 
layer, and A = d 2 /dx 2 + d 2 /dy 2 + d 2 /dz 2 the three-dimensional Laplace-operator. 
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The current-density vector jii(x, y, z, t) inside the p bulk is connected with the local 
potential <p(x, y, z, t) according to 



3r = -vr^V, 



(37) 



and its x component ]r >x at the boundary x = wr "feeds" the thyristor-like subsystem 
locally so that 



Jr. 



X I X=WR 



3t{v,z). 



(38) 



The p + -p -n-p + -n ++ sample is driven by an ideal voltage source Vg via an external load 
resistor R L . This yields the following load line equation: 



V A (t) = V S -R L [ j T (v, z, t) dydz. 



(39) 



The sets of equations (|) and <g), (|TT|), (|T|), (H) - @ present a self- 

consistent non- stationary 1 x 2-dimensional description of the system under consideration 
for the low and moderate current range, respectively. 



IV. HOMOGENEOUS STATIONARY STATES AND CORRESPONDING 
CURRENT-DENSITY VS. VOLTAGE CHARACTERISTICS 



To describe the homogeneous stationary states Eqs. (||), (|TTD, ([HI), ([Hp, and (0) for 
the low current range and Eqs. flll]), (ID, (H), ©, and (||) for the moderate current 
range with all time and space derivatives set to zero have to be solved: dQ/dt = 0, Vj_Q = 
0, Vj_y9i = 0, V±ip 2 = 0, V±ip = 0, V±V C = 0, V ±V T = 0. This leads to j T = h = jc = 
ji = j = const. 



A. Low current range 



The current gain ot\ of the n-p + -n ++ transistor can be written as 



Qti = ai(v?i) = — 
J 



j lns exp (qtfi/kT) 



(jina + hps) exp (qtpj kT) + j 1Rs exp (qcp 1 /2kT) ' 



(40) 



Combining the charge balance equation (|]) for the stationary homogeneous state with j 2n 
jcn = otij and Eq. (g^) yields 



Q = = T*jinsexp (qipi/kT), with 



= r 



(41) 



According to Eq. ( |TUD the variable Q is proportional to jc p , which is equal to (1 — a\)j; 
together with Eq. (Rj) and ~ Vc one obtains: 



Q = t' 



1/2" 



V r 



pth , 



«l) j- 



(42) 
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Combining Eqs. ( [41]) and (f42|) leads to the following parametrically determined j(Vr) 
characteristic of the thyristor-like V-n-p + -n ++ subsystem for the low current range: 



V T = V 



pth 



1 - 



(43) 



where the voltage drop (pi across the p + -n ++ emitter and the current density j are coupled 
by 



j = (jins + hps) exp (qipi/kT) + j XRs exp (qip 1 /2kT). 



(44) 



The differential resistance of the thyristor-like subsystem for the homogeneous state under 
low current conditions is then given by: 



dVp dVr da.\ dipi 
dj da i d<fi dj 



(45) 



dtpi 
dj 



dj 



x 



B. Moderate current range 



Using the charge balance equation (§4|) and jin = j — jcp = oi\ j we obtain the following 
relation between Q and j: 



with 



Q = T *9oJ, 



(46) 



6 + 1 



On the other hand, for the homogeneous case the variable Q depends on j according to the 
relation 



Q = Qr 



1/2 



V, 



pth , 



1/2 



1/2 



(47) 



which follows from Eqs. @, (|26|) together with j Cp = (1 — ai)j and Vr ~ Vc- 
Choosing 

. _ qN D w%,w 20 (l - a>i) 



WpTlgl 



(4J 



and equalizing Eqs. (|46|) and ( ^7|) yields together with Eq. ( ^7|) the following expression for 
the j(Vr)-characteristic of the thyristor-like subsystem for the case of moderate currents: 



J = Jr 



V, 



T 



1/21 



pth / 



(49) 
(50) 
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where j m and V pt h are given by Eqs. ( [48|) and (|9]). 

The specific differential resistance of the thyristor-like subsystem in the homogeneous 
state for this case is: 



with r. 



in 



2V r . 



pth 



Ji 



C. The global current-density vs. voltage characteristic 

The global i(V^) -characteristic of the complete p + -p~-n-p + -n ++ structure ensues from 

V A (j) =V t + V r , (52) 

where Vr = r R j denotes the voltage drop across the resistive layer, the resistance of which 
is determined by r R = wr/ \q^ p Nj). 

D. Calculated current-density vs. voltage characteristics 

By properly adjusting the design parameters of the n ++ layer it is possible to realize 
the appearance of a negative differential resistance for both the low and moderate current 
case. In the calculations presented below the saturation current-densities j\ ns and ji ps have 
been used to adjust the two cases. Figure |3](a) shows stationary current- density vs. voltage 
characteristics for the low current case. Besides the characteristic j(Vr) of the thyristor-like 
subsystem and the characteristic of the total structure j (Va), which are both calculated for 
the set of standard parameters, two further j(V^)-characteristics for reduced values of the 
conductivity of the p~ bulk are plotted. For standard parameter values the bulk resistance 
of the p~ bulk is not sufficiently large to compensate the negative differential resistance 
of the thyristor-like subsystem at low currents. For this reason the two curves j(Vr) and 
j'(Va) are nearly identical and cannot be distinguished in the diagram. In fact, the current- 
density jNDR, at which the negative differential resistance transforms to a positive one, are 
nearly equal and approximately 1 mA/cm 2 for the two cases. Jndr can be reduced by a 
considerable amount only, if the resistance of the p~ bulk is enlarged at least by a factor of 
ten, as illustrated by the two other curves in Fig. |3|(a) . 

Figure ^|(b) shows current- density vs. voltage characteristics for the case of moderate 
currents and the standard value of the conductivity of the p~ bulk. The current-density 
jNDR, at which the differential resistance changes from negative to positive values, is 0.45 
A/cm 2 for the total structure, i.e. considerably lower than the corresponding value (^ 1 
A/cm 2 ) for the thyristor-like subsystem. 

The blocking voltage is approximately equal to the punch-through voltage V p th for both 
the low and moderate current case. This is because the emitter efficiencies are very small 
at low currents, so that a hole injection current from the anode emitter, that is sufficiently 
large to initiate the regenerative process, becomes possible only, when the effective width 
W2{Vc) of the n base is approximately zero, i.e. under punch-through condition. 
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V. TRANSVERSAL STABILITY ANALYSIS 



In this section we analyze the stability of stationary homogeneous states along the 
transversal dimensions with respect to small 1 x 2-dimensional fluctuations of the space- 
dependent variables Y. The standard approach for the linearization procedure implies to 
make the following Ansatz for all variables in the vicinity of a stationary state Y sta t 

Y(x, y, z, t) = Y stat + 5Y(x, y, z) exp ((t). (53) 

The small variations 5Y(x, y, z) introduced here have to satisfy the same boundary con- 
ditions as the variables Y(x, y,z, t) [see Eqs. (|30|)-(|3"3|) and (p5|), (|36|)] . 



A. Fluctuations of the potential in the resistive layer 



The distribution of the potential in the resistive layer, which is essentially three- 
dimensional, can be described as solution of the Laplace equation (|34]) with the boundary 
conditions fl3"5] ) and fl3T)|). A general solution of this equation is a superposition of a uniform 
excitation mode and various /c-modes describing spatially periodic potential perturbations: 



with k„ 



(p(x,y,z,t) = 
= mir/Ly, k z 



-rr V A ~ V T 

V a x 

= n7r/L z , k 2 = 



Bk{t) sinh {kx) cos (k y y) cos (k z z), 



(54) 



k 2 + k 2 , and m,n denoting integers. Figure |] il- 
lustrates the potential fluctuations 8ip(x,y,z) as well as the potential <p(x,y,z) in a re- 
duced two-dimensional space. For a given potential fluctuation 5(p(x = WR,y,z,t) = 
5ip m cos {k y y) cos ik z z) exp ((t) at the boundary wr between the resistive layer and the 
thyristor-like subsystem the fc-dependence of the potential variations 5<p(x, y, z, t) in the 
resistive layer is given by: 

sinh (kx) 



8ip(x,y, z,t) = 5<f T . 



cos (k y y) cos (k z z) exp ((t). 



(55) 



sinh (kwji) 

The x component jr >x of the current density Jr{x, y, z, t) at the border between the 
thyristor-like subsystem and the resistive layer is equal to the local anode current-density jr 
of the thyristor-like subsystem [see Eqs. (|37|), (|38|)1, so that the following relation is valid: 



JR,. 



JT = -<7Jl a~ 
ox 



x=w R 



_ d(p 



(56) 



x=w a 



From Eqs. (|55D and (|5q ) one obtains the following variational derivative by which fluctuctions 
of 8jt and 5Vr are coupled for each wave-number k: 

5V T 



$3t 



with 



R 



wr tanh (kvjR) tanh (kwR) 

— — Ti 



(57) 



qfi p N} 



kw 



R 



R- 



kw 



R 



The quantity r R can be considered as an effective specific resistance of the bulk material with 
respect to a potential fluctuation with wave-number k. It differs from the specific resistance 
r R — w r/ {q^p^a ) °f tlie p~ bulk because of three-dimensional deformations of the potential 
relief which arise due to a non-uniform current flow in three dimensions. 
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B. Fluctuations in the thyristor-like subsystem 

Now we apply the linear stability analysis to the charge dynamic equation making a 
perturbation Ansatz for all the variables Y(y, z, t) subjected to two-dimensional transversal 
fluctuations, 

Y(y, z, t) = Y stat + 5Y m cos(k y y) cos(k z z) exp(Cfc), (58) 
and then study the stability for the both cases mentioned above. 



1. Low current range 



For the low current range the following set of variational relations can be derived from 
Eqs. d)-© and 



8j2n 

SQ 
where 



C + - + k z D n ) 8Q, 

r* 



a k c + 



dQ 



l-a k c 



R 



dVr 



r k R 5j T + 



dQ 



stat 



l Cj 



1 + 



djcp 
'l_ 1_ 



5j T + 



a k c k 2 D n -(l-a h c )^-^— 

H < 2 ^Cstat 



5Q, 



stat 



R 



6j T + (1 - a k c ) 



k 2 D n + 



kT 1 



Q r 2 Qstat 



5Q, 



dQ 



dV T 



stat 



a c — oi\ 1 
1 — oti r>p 



dQ 



dj 



c P 



aci(ipi 



stat 



1 - ai(<pi)' 



Qstat = T*a>ij T , 



(59) 
(60) 
(61) 
(62) 

(63) 



and the differential resistance of the thyristor-like subsystem for the homogeneous state is 
given by: 



1 



7* ai 
t' 1 — ai 



(64) 



The coefficient a c is defined as ratio of the fluctuations of the electron current-density 
5jcn an d the total current density Sjc at the collector, which in turn can be written as: 





dji n / d<pi 




dji/dipi 




_ dji n /dipi 


Sji 


dji/dipi 



' 1 / ' 



era 

In 



dji/dipi \1 + rf m /r k J l+rf m /r k 



for k 0, and 



(65) 



' In 



The differential resistances rf m and r^™ of the p + -n ++ emitter junction are defined as 

dji 



„em\ — 1 



,em\ — 1 



In 



difi 
djln 
difi ' 



(66) 
(67) 
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The fc-dependent quantities r k , r k , and r\ can be considered as effective resistances of the 
p + base, the n base, and the V layer, respectively, for a current-density fluctuation with 
wave-number k. 



(r k ) 1 = qii p N\ Wl k 2 , 



(r k 2 (h 



qH n N D w 2 {V T ) 



k . 



(68) 
(69) 

(70) 



stat 



Using these dependencies we obtain from Eqs. (|59"D-(|6"2) the final dispersion relation for 
the increment ( as a sum of four components: 



C{k,j T ) = Co + G 
where the partial components are: 



C2 + C: 



Co 

Ci 

Cv 
Cv 

Cv 



a£ — a\ 
T*K{h) 

ac r i m 
rJC{k) 

C2 ~l~ C2 > 



X 



„/>- ■ 



(71) 

(72) 

(73) 
(74) 



r*/C(&) 



x 



'c 



r*/C(/c)«i 



-k 2 D r 



a 1 



K{k) 



r fc> 



1 + ^ 



Ql 1 + 



K a c - a h 



kT 
x x 

QJT 



R 



,k ■ 



(75) 



'v. 



{a k c - ai )r k R 
T*JC(k) 



x 



The function K.(k) is given by the expression: 



)C(k) =ai(l-a k c ) 



r k 



K a c - a lt 



R 



(76) 



Note that all partial components depend on the current- density dependent coeffficients a°, 
a k and at\. Besides, the contribution of the terms d, (2, and £p is controlled by the effective 
resistances of the p + base, the n base, and the V layer, respectively. For (-p also the diffusion 
properties of the plasma charge play an important role. 



2. Moderate current range 

For the moderate current range a corresponding set of variational relations can be ob- 
tained by applying the linearization procedure to Eqs. (p0| ) -(|29| ) with taking into account 
the relations Eqs. fliDD and (|57|): 
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Sj 2n = U + - + k 2 bD p ) 5Q 



6 + 1 



1 + Uh 



'V. 



Sj 



2n 



SQ 



<&+[-% 



l-a k c 



'v 



h 

4 



dQ 



dV T 
Sjc P = (1 - 4 
where 



r k R 5j T + 



dQ 



stat 



djcp 



Sj T + 
Sjc P , 



k 2 a k c (b - 1)D P - (1 - a k c ) 



kT 1 

£[Qstat ^2 



SQ, 



stat 





5j T + 


\+(- + -)r k ' 


\r k Tp) R 





(l-a k c )(b-l)k 2 D p + (l-a 



c, 



kT 1 

qQstat r\ 



dQ 



C \ _|_ j-em I yk ' 



dV T 
dQ 



stat 2r^ 



djcp 



stat 



2(1-0;/ 



(77) 
(78) 
(79) 
SQ, (80) 

(81) 
(82) 
(83) 



The parameter can again be considered as an effective resistance of the plasma layer with 
respect to current-density and potential fluctuations with a wave-number k. However, in the 
case of moderate currents is conductivity modulated: 



rp(Qstat) = ((b + l)fi p Qstat(jT)k 2 ) 1 , 

Qstatiir) = T*g j T - 

The final dispersion relation for the increment ( results from Eqs. (|77|)-(|80 

CiKh) =Co + Cl + C2 + G>, 

with 

rJC(k) { r°J' 



Co 

Ci 



(84) 
(85) 

(86) 
(87) 



'2(1 -ai) 
9o 



a c r i™ 



X 



C2 = C2 + C2 1 

( 2{l-a x ) \ (l-a k c )r k R 1 

Q ~ { 9 o + ) rX{k) X rl 



.k ' 



39) 



k ■ 

2 



C" 

S2 



9k\l + 



(1 - a£)4 + K{k) 



l-a k c kT 1 
T*g lC(k)qj T rf 



Cr — Cv + Cp> 



-fc 2 L> p i 6(1 - a k c ) + a k J + 



(ft-l)(l-q*) 
/C(Jfe) 



(90) 



k \ R 



2g k (l - ax) - g (l - a 



x 



ng JC{k) " r%' 
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and 

9k = a k c ~ 



b 

b + r 



(91) 



(fytk /yi fa \ ty>fa 

l + lf + lf -(1-ax)-? >0. (92) 
'2 '-p/ ' T 

Thereby the differential resistance of the thyristor-like subsystem is determined by 
Eq. (P. 

For both the low and the moderate current case the increment ((h) is proportional to the 
sum of the differential resistance of the thyristor-like subsystem, r^, and that of the resistive 
layer, r^, in the limit of perturbations with small wave-numbers: 

C(*-0)oc-(r£ + r°). (93) 

This obviously corresponds to a one-dimensional loading of the thyristor-like subsystem by 
the resistive layer. From Eq. ( |93| ) follows the well-known criterium for stability with respect 
to homogeneous fluctuations: 

r< k + r T> 0- (94) 



C. Dispersion relations and discussion 

Figure [| shows dispersion relations ((h) for the low current case. The calculations have 
been performed with two reduced values of the p + base conductivity al in order to illustrate 
two qualitatively different instabilities that can occur in the semiconductor structure. If the 
p + base conductivity is sufficiently low, dispersion relations like the curves denoted by (2) 
can occur. They are characterized by a pronounced maximum at a critical wave-number 
that is definitely larger than zero. Curve (2a) represents the case that the increment ( of 
this critical wave-number is zero, while all other increments are still smaller than zero, i.e., 
the uniform state becomes destabilized with respect to current-density fluctuations with the 
critical wave-number. When the sample current is decreased, the curve is shifted up [curve 
(2b)], and the increments of a certain band of wave-numbers are larger than zero. This 
behavior can be interpreted as a Turing-like instability^. 

When the transversal coupling inside the thyristor-like subsystem is strengthened, e. g. 
by reducing the transversal resistance of the p + base, the maxima of ((k) shift to smaller k 
values leading to dispersion relations like those denoted by (1). They are characterized by 
a monotonous decrease of ((k) with growing k, which can be considered as an instability of 
Ridley's typeEl. In such a case, the uniform state is typically destabilized via a saddle node 
bifurcation by fluctuations with the shortest possible wave-number 2ir/(Ly + L^), if uniform 
fluctuations are suppressed by a sufficiently large load resistor in the external circuit. Similar 
to the Turing-like case, the curve is shifted upwards, when the sample current is decreased. 
Note that the parameter j is the same for the curves (la) and (2b) and has been adjusted 
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such that the differential resistance of the global current-density vs. voltage characteristic 
vanishes, i. e., + = 0. 

Dispersion relations ((k) for the moderate current case are shown in Fig. |6]. Similar to 
the low current case the homogeneous state is destabilized at a critical current level either by 
an instability of Ridley's (curves 1) or Turing's type (curves 2). The bifurcation type can be 
controlled by a proper design. The most important parameters are the conductivities of the 
resistive layer, the p + and the n base. The former essentially governs the transversal spread- 
ing of the potential drop across the resistive layer and tends to damp out current-density 
fluctuations. The p + and n base conductivity control the transversal hole and electron cur- 
rent densities in the two base layers. As these base currents in turn regulate the injection 
of electrons and holes from the n + emitter and the plasma layer V, respectively, they play 
an important role concerning the transversal spreading of a region with enhanced current 
density. The competition between the damping and the activating process, is the essential 
physical mechanism allowing - under properly chosen design parameters - the destabilization 
of a homogeneous current-density distribution by spatially periodic current-density fluctua- 
tions strongly reminding of the generic mechanism proposed by Turing. 

VI. SUMMARY 

In conclusion we point out that the suggested theory, which is based on a subdivision of 
the semiconductor structure into an active and a passive subsystem, supplies a self-consistent 
quantitative description of the nonlinear mechanisms, which control the longitudinal current 
flow, and the transversal current-density instabilities in multilayer semiconductor systems. 
The nonlinear longitudinal electrical properties in the active layer are based on a thyristor- 
like regenerative mechanism, while the passive subsystem acts as a simple resistive layer. 
Depending on the system parameters two different types of regenerative mechanisms can 
be distinguished leading to specific model equations. A stability analysis of the derived 
equations reveals that in each case two basic bifurcation types can be expected when the 
total sample current is varied. The bifurcations are caused by instability mechanisms that 
are closely related to those studied by Turinga and RidlejM. It turned out that in particular 
the conductivity ratios of the p + base, the plasma layer, and the resistive layer are important 
parameters to control the instability mechanism. Thus, adjusting the respective transversal 
conductivities, e.g., by irradiating the sample with protons or ions, using compensated 
substrate material, or adapting the sample geometry, it is possible to design a multilayer 
system with the aspired (in) stability features. Finally, we emphasize that the suggested 
approach presents a basis for further numerical analysis of pattern formation processes. 
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TABLE I. Physical constants and material parameters 



Symbol 


Symbol meaning 


Value 


b 


b = D n /D p 


2.8 


D n 


diffusion coefficient of electrons 


D n = kTfi n /q 


D p 


diffusion coefficient of holes 


Dp = kTfip/q 


k 


Boltzmann constant 


8.62 • 10 S V/K 


Hi 


intrinsic carrier concentration 


1.45 • 10 10 cm" 3 


q 


elementary charge 


1.6 • 10~ 19 As 




vaccuum permittivity 


8.854 • 10~ 14 F/cm 


e 


dielectric constant of Si 


11.7 




electron mobility 


bfip 




hole mobility 


480 cm 2 /Vs 


TABLE II. Set of design parameters 


Symbol 


Symbol meaning 


Value 


Jlns 


saturation value of the electron injection current-density at 


2 • 10- 11 Acm^ 2 




the p + -n ++ junction 




Jlps 


saturation value of the hole injection current-density at the 


1 • 10~ 13 Acm- 2 




p + -n ++ junction 




jlRs 


saturation value of the recombination-generation current- 


1 • 10~ 8 Acm" 2 




density at the p + -n ++ junction 




N A 


acceptor concentration of the p + emitter 


1 • 10 18 cm" 3 


N A 


acceptor concentration of the p~ bulk 


1 • 10 13 cm" 3 


"I 


acceptor concentration of the p + base 


5 • 10 17 cm~ 3 


N D 


donor concentration of the n base 


1.2 • 10 14 cm" 3 




donor concentration of the n ++ emitter 


1 • 10 20 cm" 3 


W\ 


width of the p + base 


2 • 10~ 4 cm 


W 2 Q 


width of the n base for Vq = 


1.5 • 10~ 3 cm 


W-p 


width of the plasma layer 


5 • 10~ 3 cm 


WR 


width of the resistive layer 


1 cm 


(?\ 


conductivity of the p + base 


38.4(Ocm)- 1 




conductivity of the n base 


2.58 • l®- 2 {VLcmy l 


a-p 


conductivity of the V layer 




OR 


conductivity of the p~ bulk 


7.68 • 10- 4 (ficm)- 1 


r 


lifetime in the V layer for low currents 


1/xs 


Th 


lifetime in the V layer for moderate currents 


5/Us 
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FIG. 1. Geometry of the sample (a) and doping profile along the vertical (anode-to-cathode) 
direction (b). 
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FIG. 2. Non-equilibrium plasma-field stratification along the vertical direction (a) and possible 
contours of small current-density fluctuations which are periodical along the transversal directions 
(b). 



1.0 

0.8 

la o.6 

& 0.4 
0.2 
0.0 

5 10 15 20 25 5 10 15 20 25 

voltage (V) voltage (V) 

FIG. 3. Current-density vs. voltage characteristics for the low (a) and moderate current range 
(b). Va and Vt denote the bias of the total system and the thyristor-like subsystem, respectively. 
Parameter values for (a) are indicated in Table O, modified parameters for (b) are wr = 2 • 10~ 2 
cm and ai = ji ns /{jins + hps) = 0.821. 
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FIG. 5. Dispersion relations C(k) for the low current case revealing the instability of the uniform 
state with respect to large wave-length fluctuations, k — > 0, (1) and with respect to a fluctuation 
with a critical wave- number k > (2); a\ SS = 0.04(71,5 • 10~ 4 ai, for curves (1) and (2), and 
j = 1.01,0.671,1.11,1.01 mA/cm 2 for curves la, lb, 2a, 2b, respectively. Other parameters as in 
Table 0. 
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FIG. 6. Dispersion relations <^(k) for the moderate current case revealing the instability of the 
uniform state with respect to large wave-length fluctuations, k — > 0, (1) and with respect to a 
fluctuation with a critical wave- number k > (2); af f = 0.25cti,0.05 • lO" 4 ^, for curves (1) and 
(2), and j = 394,350,407,394 mA/cm 2 for curves la, lb, 2a, 2b, respectively. Other parameters 
as in Table [D] besides ./V^ = 5 • 10 13 cm~ 3 and wr = 0.1 cm. 
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